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CONVERSION  FACTORS,  NON-SI  TO  SI  (METRIC)  UNITS  OF  MEASUREMENT 


Non-SI  units  of  measurement 
(metric)  units  as  follows: 

_ Multiply _ 

cubic  feet  per  second 
feet 

knots  (international) 
miles  (US  nautical) 
miles  (US  statute) 
square  miles  (US  statute) 


used  in  this  report  can 

_ g y. 

0.02831685 

0.30*18 

0.5144444 

1.852 

1.609347 

2.589998 


be  converted  to  SI 

_ To  Obtain 

cubic  metres  per  second 
metres 

metres  per  second 
kilometres 
kilometres 
square  kilometres 


PORT  OF  MIAMI  NUMERICAL  MODEL  STUDY 


PART  I :  INTRODUCTION 

Background 

1.  The  Port  of  Miami,  Florida,  is  located  in  the  northern  part  of  Bis- 
cayne  Bay  and  along  both  banks  of  the  lower  part  of  the  Miami  River  (Fig¬ 
ure  1).  The  major  deepwater  navigation  entrance  to  the  port  is  via  a  dredged 
channel  extending  several  miles  from  the  ocean,  through  Government  Cut,  and 
across  the  bay  to  the  turning  basin  at  the  municipal  terminal.  Because  of 
problems  with  strong  cross  currents  at  the  ocean  entrance  and  inner  harbor 
channels  to  Miami  Harbor,  the  US  Army  Engineer  District,  Jacksonville  (CESAJ), 
is  evaluating  Jetty  modification  plans  at  the  entrance  channel.  To  better 
quantify  the  problems  and  to  evaluate  remedies,  the  US  Army  Engineer  Waterways 
Experiment  Station's  (CEWES's)  Coastal  Engineering  Research  Center  (CERC),  was 
requested  to  conduct  a  tidal  circulation  numerical  model  study. 

2.  CERC  applied  the  WES  Implicit  Flooding  Model  (WIFM)  (Butler,  in 
preparation),  a  two-dimensional  (2-D)  vertically  integrated  numerical  model, 
along  with  a  companion  prototype  data  collection  study  in  the  Miami  area.  The 
field  data  required  for  numerical  model  calibration  and  verification  were  col¬ 
lected  by  the  National  Oceanic  and  Atmospheric  Administration  (NOAA)  and  CERC. 
All  numerical  computations  were  performed  on  the  Control  Data  Corporation's 
CYBER  205. 

3.  The  numerical  modeling  approach  consisted  of  applying  WIFM  in  a  two- 
phase  modeling  study.  The  first  phase  involved  a  global  model  which  included 
areas  of  the  inner  harbor  and  a  substantial  open  ocean  area  with  a  200-  to 
300-ft*  minimum  grid  cell  dimension.  The  horizontal  grid  spacing  was  vari¬ 
able,  and  the  finer  resolution  was  concentrated  at  the  ocean  entrance  and  port 
facilities.  This  model  was  calibrated  against  neap  tidal  events  and  was  veri¬ 
fied  against  spring  tidal  events  using  measured  data. 

4.  The  second  phase  of  the  modeling  effort  involved  the  development  of 
a  fine  resolution  grid  (minimum  grid  dimension  of  100  ft)  in  the  entrance 

*  A  table  of  factors  for  converting  non-SI  units  of  measurement  to  SI 
(metric)  units  is  presented  on  page  3. 


channel  and  port  facility  areas.  Boundary  conditions  for  this  grid  were 
supplied  by  the  global  model.  This  model  was  calibrated  and  verified  for 
tidal  events  using  observed  data.  The  refined  numerical  model  was  necessary 
to  permit  accurate  determination  of  the  impact  of  small-scale  changes  to  the 
port  and  entrance  channel  system  on  the  coastal  current  regime  and  port 
hydrodynamics. 

Scope 

5.  This  report  describes  the  use  of  the  global  and  refined  numerical 
models  to  evaluate  the  tidal  current  regime  at  the  entrance  to  the  Port  of 
Miami  and  to  assess  the  impact  of  structural  (Jetty  extensions)  and  nonstruc- 
tural  modification  plans  for  reducing  strong  cross-current  effects  on  navi¬ 
gation.  Hydrodynamic  simulations  using  the  refined  grid  model  were  made  to 


provide  input  data  into  the  vessel  simulation  study  conducted  by  the  US 
Department  of  Transportation  Maritime  Administration's  Computer-Aided  Opera¬ 
tion  and  Research  Facility  (CAROF).  This  study  was  sponsored  by  CESAJ  for  the 
Port  of  Miami. 


PART  II:  CHARACTERISTICS  OF  THE  STUDY  AREA 

Study  Area  Description 

6.  The  study  area,  Port  of  Miami,  Florida,  is  located  in  the  northern 
part  of  Biscayne  Bay  and  along  both  banks  of  the  lower  6-mile  reach  of  the 
Miami  River  (Figure  1)  which  empties  into  the  west  side  of  the  bay.  The 
northern  part  of  the  Bay  is  separated  from  the  Atlantic  Ocean  by  Miami  Beach 
Peninsula  and  the  southern  part  by  Fisher  Island,  Virginia  Key,  and  Key 
Biscayne.  Average  water  depths  in  Biscayne  Bay  range  from  5  to  10  ft  at  mean 
low  water  (mlw) . 

7.  The  main  mechanism  responsible  for  flow  in  and  out  of  the  bay  is  the 
semidiurnal  astronomical  tide.  The  mean  neap  tide  is  about  0.8  ft  (above  the 
National  Geodetic  Vertical  Datum  (NGVD)),  and  the  mean  spring  tide  is  about 
1.3  ft  in  the  bay.  Tidal  ranges  in  the  bay  are  about  81  percent  of  the  tidal 
range  in  the  open  coast.  Mean  tide  variation  between  the  harbor  entrance  and 
th*>  bay  is  about  0.2  ft.  Strong  easterly  winds  raise  the  water  level  about 
0.7  ft  at  the  entrance  to  the  port  and  0.5  ft  in  the  bay,  while  strong  west¬ 
erly  winds  lower  the  water  level  about  0.5  ft  at  the  entrance  and  about 

0.25  ft  in  the  bay  (Headquarters,  US  Army  Corps  of  Engineers  1982). 

8.  The  major  deepwater  navigation  entrance  to  the  port  is  via  a  dredged 
channel  extending  several  miles  from  the  ocean,  through  Government  Cut,  and 
across  the  bay  to  the  turning  basin  between  Fisher  and  Lumus  islands.  Gov¬ 
ernment  Cut  is  flanked  on  either  side  by  two  stone  Jetties  (Figure  1),  one 
extending  from  the  southern  end  of  Miami  Beach  and  the  other  from  the  north¬ 
easterly  side  of  Fisher  Island.  On  the  west  side  of  this  island  are  two  man¬ 
made  islands,  triangular-shaped  Lumus  Island  and  rectangular-shaped  Dodge 
Island,  the  main  function  of  which  is  to  handle  cargoes. 

Ship  Maneuverability  Problems 

9.  Large  ships  sailing  into  the  Port  of  Miami  encounter  problems  with 
cross  currents  at  two  distinct  locations  along  the  dredged  channel.  These  two 
sites  were  identified  by  the  ships'  pilots  as  locations  A  and  C  (Figure  1). 

The  problem  at  location  A  occurs  when  northeasterly  winds  blow  over  area  A. 
These  winds  drive  buoy  2  into  the  main  channel,  as  shown  in  Figure  1.  Ships 
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sailing  toward  the  harbor  are  forced  to  maintain  a  crab  angle  about  4°  S 
(according  to  pilots).  Visible  evidence  on  Gulf  Stream  currents  also  has  been 
reported  on  the  southern  side  of  the  dredged  channel.  The  main  edge  and  width 
of  the  Gulf  Stream  are  unknown.  There  appears  to  be  no  literature  on  Gulf 
Stream  migrations  in  the  vicinity  of  the  Port  of  Miami;  however,  pilot  reports 
indicate  that  the  Gulf  Stream  edge  migrates  westward  as  far  as  channel  Buoy  7A 
along  the  channel. 

10.  Vessels  leaving  the  dogleg  turn  (position  B)  enter  Government  Cut 
along  the  center  line  of  the  channel.  During  their  sail  they  must  maintain  a 
course  toward  the  south  jetty  to  avoid  collision  with  Buoy  8  (Figure  1)  which 
penetrates  into  the  main  channel  when  northeasterly  winds  blow  over  the  area. 
As  ships  come  abreast  of  the  north  jetty,  they  must  maintain  a  crab  angle  to 
avoid  collision  with  Buoy  9  (position  C).  This  near-collision  is  caused  by 
northeasterly  winds  pushing  the  stern  of  the  ship  toward  the  south  jetty, 
thereby  causing  the  bow  to  move  toward  the  north  jetty.  Northeast  wind 

effe ‘ts  are  partially  blocked  by  the  north  Jetty  as  the  ship's  stern  passes 
the  tip  of  the  north  jetty.  At  this  time  the  vessel  is  pushed  toward  the 
north  jetty  (position  D)  by  strong  cross  currents  generated  at  the  tip  of  the 
south  Jetty. 

11.  As  southwesterly  winds  and  cross  currents  intensify  the  movement  of 
Buoy  9  into  Government  Cut  Channel,  vessels  must  follow  a  course  closer  to  the 
north  Jetty  to  avoid  colliding  with  the  buoy.  Under  normal  wind  conditions 
ship  collision  is  not  a  problem;  however,  if  the  wind  and  current  forces  on 
the  ship  are  increased,  the  vessel  can  be  pushed  closer  to  the  north  jetty  and 
possibly  grounded.  Pilot  reports  indicate  grounding  of  ships  at  the  north 
jetty  under  gusty  southwesterly  winds. 

12.  During  flood  tide  with  northeasterly  winds,  large  vessels  avoid 
meeting  each  other  in  the  vicinity  of  the  Jetties  for  the  reasons  described 
above.  This  restriction  prevents  timely  arrival  and/or  departure  of  ships  and 
endangers  ship  safety.  A  similar  but  less  serious  situation  occurs  when  ves¬ 
sels  leave  the  port  during  ebb  tides. 

13.  In  its  application  of  WIFM  to  these  problems,  CERC  will  (a)  eval¬ 
uate  tidal  current  regime  at  the  entrance  to  the  Port  of  Miami,  (b)  assess  the 
impact  of  proposed  jetty  extensions  on  reducing  strong  cross-current  effects 
on  navigation,  and  (c)  evaluate  the  impact  of  proposed  port  facility  changes 
on  tidal  circulation. 
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PART  III:  COMPUTATIONAL  METHODS 


Governing  Equations 


14.  The  governing  equations  WIFM  uses  for  hydrodynamic  calculations  are 


as  follows: 


Continuity  Equation 


ft  *  if  (ud)  +  if  (vd)  =  R 


Momentum  Equation  ( x-direction ) 


3u  3u  3U  c.  3  ,  ,  gu  ,  2  2. 

it  *  u  ii  *  7 1?  '  fv  * 8  il (n  *  \}  *  ^  <u  *  v  > 


-  e(^  ,  ifu),  F  .  o 

lax2  ay2/  * 


Momentum  Equation  (y-direction) 


3v  3 v  3v  „  3  ,  ,  g v  ,  2  2, 

3t  *  U  ii  *  V  5J  *  fu  *  8  i?  (,>  '  "a’  *  ^  (u  *  v  > 


(>2V  3^^  p  0 

e  l  — r  +  — »  /  +  F  =0 

\3x2  3y2/  y 


in  which  n*  ,  u  ,  and  v  are  the  dependent  variables  representing  water 
surface  elevation  above  datum  and  vertically  integrated  velocities  in 


*  For  convenience,  symbols  and  abbreviations  are  listed  in  the  Notation 
(Appendix  A). 


the  x-  and  y-directions,  respectively.  Independent  variables  in  the  above 

equations  are: 

d  =  n  -  h,  total  water  depth 
h  =  bed  elevation  above  datum 
t  =  time 

f  =  Coriolis  parameter 

g  =  acceleration  due  to  gravity 

n  =  hydrostatic  water  elevation  due  to  atmospheric  pressure 
differences 

C  =  Chezy  coefficient 

e  =  eddy  viscosity  coefficient 

FX,FV  =  external  forces  in  the  x-  and  y-directions,  respectively 
*  (  i.e. ,  wind  stress) 


Numerical  Formulation 

15.  WIFM  uses  the  alternating-direction-implicit  (ADI)  scheme  to  solve 
Equations  1-3.  However,  because  of  the  inclusion  of  the  advective  terms  in 
these  equations,  the  ADI  scheme  encountered  stability  problems.  To  minimize 
these  problems,  WIFM  uses  a  centered  stabilizing-correction  (SC)  scheme,  that 
is  accurate  to  second  order  in  space  and  time.  Boundary  conditions  can  be 
formulated  to  the  same  order  of  accuracy.  A  brief  description  of  this  tech¬ 
nique  is  given  below.  Butler  (in  preparation)  provides  additional  details  on 
the  SC  scheme. 

16.  Equations  1,  2,  and  3  can  be  written  in  matrix  form  as 


U.  +  AU  +  BU  =0 
t  x  y 

where 


o  d  o 

A  -  goo, 
0  0  0 


o  o  d 


goo 


The  SC  scheme  for  solving  Equation  4  is 


(1  +  x..)  U»  =  (1  -  x._  -  2X..)  U 


(1  +  X  )  U  =  U*  +  x  U  (6) 

y  y 

where 

X  =  ^  7^  A6  ,  2-D  difference  operator 
x  2  Ax  x  ’  K 

At  =  time-step 

Ax  =  length  of  computational  cell  in  x-direction 
6X  =  centered  difference  operator 

U*  =  Value  of  u  at  an  intermediate  time-step  level  (k  +  1  time  level) 
k  =  integer  time-step  counter 
*  =  intermediate  time-step  level 

X  =  ^  ^  B6  ,  2-D  difference  operator 
y  y  y 

Ay  =  length  of  a  computational  cell  in  y-direction 
6^  =  centered  difference  operator 

17.  The  SC  scheme  consists  of  two  steps:  approximating  the  grid  in 
the  x-direction  and  sweeping  the  grid  in  the  y-direction.  Completing  both 
sweeps  constitutes  a  full  time-step  and  marches  the  solution  from  the  tim 
level  to  the  k  +  1  time  level.  The  forms  of  the  continuity  and  momentum 
equations  employed  in  the  multioperational  hydrodynamic  scheme  are  given  by 
X-sweep 


1  ,  „  k-K  1 

<n  -  n  >  *  24 


—7—  &  (u^d+u'^d)  +7—6  (v<-^d)  = 

2Ax  x  Ay  y 


IE  <u*  -  u  >  *  0 


24t  <»*  -  '"'l  'Zf  6y  -  0 


Y-sweer 


1  ,  <-1 
2At  (n 


K+  1  , 

u  =  u* 


"*>  *  2iy  5,  (v<‘lcl  -  v<’’d) 


1  ,  k+1  e.  ,  <+ 1  <- 1 . 

IE  <v  -  v  >  *  2??  5y  <n  *  "  >  =  0 


18.  Notably,  v*  in  Equation  9  represents  a  functional  value  computed 
til 

for  the  k  -  1  time-level.  If  the  value  of  v*  from  Equation  9  and  the 
value  of  u*  from  Equation  11  are  substituted  into  Equations  7  and  8,  the 
following  simplified  equations  are  obtained: 

X-sweep 


^77  (n*  -  n*-1)  +57-6  (uK+1d  +  uK_1d)  +  6  (v<_1d)  =  0  (13) 

2At  2Ax  x  Ay  y 


1  ,  ic+1  tc-1%  g  .  ,  M  k-1\  « 

JKt  (U  -  U  >  +  2^  6x  (ti  +  n  >  =  0 


(14) 


Y-sweep 


I  ,  K+1  MX  1  ,  ,  <-*-1  ,  IC-1  ..  _ 

2At  71  “  n  }  +  2A7  6y  (V  d  "  V  d)  =  0 


(15) 


1  »  K+1  K-1.  g  .  /  <+1  K-1,  ft 

2fit  V  -  v  >  *  2%}  5y  (n  *  n  >  *  0 


(16) 


Equations  13  and  14  and  15  and  16  are  alternatively  solved  in  WIFM  by  applying 
these  finite  difference  equations  to  one  column  (X-sweep)  or  row  (Y-sweep), 
respectively,  of  the  numerical  grid.  Butler  (1984)  gives  the  solution  method. 

Grid  Schematization 

19.  WIFM  uses  a  stretching  transform 

x  =  a  +  b  Zc  (17) 


where 

x  =  physical  distances 
a,b,c  =  arbitrary  constants 
Z  =  computational  distance 

for  mapping  distances  along  the  coordinate  directions.  A  detailed  description 
of  the  program  MAPIT,  which  maps  a  variable  grid  in  real  space  into  a  uniform 
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grid  in  computational  space,  is  described  by  Butler  (in  preparation).  MAPIT 
maps  each  coordinate  direction  independently  and  maximizes  grid  resolution 
(finer  cells)  in  areas  of  hydrodynamic  importance  and  minimizes  the  number  of 
computational  cells  (coarser  cells)  in  areas  of  less  importance. 

Global  grid 

20.  A  global  grid  was  constructed  to  model  the  hydrodynamic  impact  of 
the  large  open  ocean.  Since  it  is  not  practical  to  make  numerical  computa¬ 
tions  in  real  space,  real  spaces  were  mapped  into  computational  spaces  by 
using  Equation  17.  This  technique  generated  a  variably  spaced  grid  which 
included  areas  of  Biscayne  Bay  north  and  south  of  Lumus  Island,  the  entrance 
channel  and  Norris  Cut,  and  approximately  40  square  miles  of  open  ocean  around 
the  entrance.  There  were  80  computational  cells  in  the  longitudinal  direction 
and  50  cells  in  the  transverse  direction.  This  grid  approximates  the  entire 
area  of  interest  with  4,000  cells.  Finer  cells  (300-ft  minimum  cell  dimen¬ 
sion)  were  provided  in  areas  of  hydrodynamic  importance  and  coarser  cells  in 
areas  of  less  importance.  Figure  2  shows  a  portion  of  the  global  grid  laid 
over  N0AA  Chart  11467  (US  Department  of  Commerce  1984). 

Refined  grid 

21.  The  second  phase  of  the  numerical  modeling  effort  consisted  of  the 
development  of  a  fine  resolution  grid  (150-ft  minimum  cell  dimension)  in  the 
entrance  channel  and  port  facility  areas.  The  finer  resolution  permits 
accurate  determination  of  the  impact  of  small-scale  changes  to  the  port  and 
entrance  channel  system  on  the  coastal  current  regime  and  port  hydrodynamics. 
This  grid,  which  includes  the  inlet  and  port  facility  areas  described  for  the 
global  grid  and  approximately  4  square  miles  around  the  entrance  to  the  port, 
is  composed  of  83  cells  in  the  longitudinal  direction  and  58  cells  in  the 
transverse  direction.  It  covers  the  area  of  interest  with  4,814  computational 
cells.  Figure  3  presents  the  refined  grid  laid  over  N0AA  Chart  11467. 

Grid  coupling 

22.  Using  the  embedded  grid  concept  to  transfer  hydrodynamic  informa¬ 
tion  from  the  global  grid  to  the  refined  grid,  WIFM  required  coupling  between 
global  grid  cells  and  other  boundary  cells  (open  ocean  boundary)  of  the  re¬ 
fined  grid.  Coupling  of  one  grid  to  another  is  done  such  that  the  cell  size 
is  the  sane  in  the  coupling  cells  of  both  gr^ds.  Leenknecht,  Earickson,  and 
Butler  (1984)  provide  additional  requirements  on  coupling  grids. 


23.  The  hydrodynamic  data  transfer  between  global  and  refined  models 


consists  of  n  ,  u  ,  and  v  at  each  coupling  cell.  In  the  present  study, 
global  data  were  saved  at  each  coupling  cell  and  obtained  for  use  in  the  re¬ 
fined  grid  by  means  of  a  utility  program  in  WIFM.  It  should  be  noted  that 
grid  coupling  is  not  totally  dynamic  in  the  sense  that  hydrodynamic  simula¬ 
tions  in  the  coupling  cells  of  both  grids  are  independent  and  not  concurrent 
Data  transfer  is  unidirectional  from  the  global  to  the  refined  model. 


PART  IV:  PROTOTYPE  DATA 


24.  Field  measurement  data  required  for  numerical  model  calibration 
and  verification  were  collected  by  NOAA  and  CERC.  Following  is  a  summary  of 
prototype  data  used  in  this  study. 

NOAA  Current  Data  and  Analysis 

25.  NOAA  conducted  an  intensive  field  data  collection  effort  in  the 
Miami  area  which  included  19  stations  located  throughout  the  study  area  (Fig¬ 
ure  4).  The  data  collection  periods  are  shown  in  Figure  5.  Current  speed  and 
direction  were  recorded  at  15-min  intervals.  These  raw  data  were  obtained 
from  NOAA,  and  a  data  base  was  constructed  from  it  for  data  analysis. 

26  A  low-pass  filtering  technique  was  performed  for  all  prototype  cur¬ 
rent  and  direction  data.  This  technique  was  necessary  to  remove  data  spikes 
caused  by  vibration  and/or  other  extraneous  behavior  of  the  instrument  and 
high  and  low  frequency  trends.  In  addition,  all  observed  data  were  analyzed 
separately  for  the  east-west  and  north-south  velocity  components  which  cor¬ 
responded  in  orientation  to  the  numerical  model  current  component.  This 
separate  analysis  allowed  for  easy  comparison  between  observed  and  predicted 
currents. 

27.  The  analysis  indicated  that  several  of  the  19  current  stations  were 
affected  by  the  instrumentation  and  calibration  errors.  Therefore,  these  sta¬ 
tions  were  discarded  from  the  analysis.  Plates  1-3  show  the  magnitude  and 
direction  of  current  recorded  at  stations  which  were  not  affected  by  the 
instrumentation  and  calibration  errors. 

CERC  Current  Data  and  Analysis 


28.  To  supplement  NOAA’s  effort,  WES  collected  current  data  from 
in  situ  current  meters  deployed  at  four  locations:  (a)  south  side  of  Gov¬ 
ernment  Cut  channel,  (b)  Norris  Cut,  and  (c)  two  locations  in  the  coastal  area 
north  of  the  ocean  entrance.  Figure  4  presents  the  locations  of  current  meter 
stations.  Table  1  shows  CERC's  current  data  collection  period  and  current 
meter  position  and  identification.  Data  collection  periods  relative  to  NOAA's 
effort  are  shown  in  Figure  5. 
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29.  Additional  current  measurements  were  taken  by  a  boat  survey  team  at 
three  ranges  in  the  entrance  channel  and  at  other  locations  throughout  the 
harbor  (Figure  4).  Some  of  the  NOAA  current  station  locations  coincide  with 
some  of  those  where  CERC  took  hand-held  current  measurements.  Overlap  of 
current  meter  locations  allowed  verification  of  the  prototype  data  accuracy. 

30.  Because  of  the  failure  of  a  number  of  NOAA's  current  meter  sta¬ 
tions,  a  considerable  amount  of  additional  data  analysis  was  required  of  both 
the  NOAA  and  CERC  data  in  order  to  develop  a  common  data  base  for  numerical 
model  calibration  and  verification.  The  analysis  and  reduction  of  all  raw 
data  were  performed  in  a  way  similar  to  the  techniques  used  for  NOAA's  current 
data  analysis.  Plates  4-8  present  the  magnitude  of  velocity  and  direction  for 
field  current  measurement  locations  76,  77,  78,  79,  and  80  (Figure  4). 

31.  Hand-held  current  measurement  data  also  were  analyzed  for  use  in 
numerical  model  verification.  Figure  6  shows  the  magnitude  of  current  mea¬ 
sured  at  selected  hand-held  current  meter  stations. 

CERC  Tide  Data  and  Analysis 


32.  CERC  carried  out  a  prototype  water  surface  elevation  data  collec¬ 
tion  program  to  obtain  field  measurements  required  for  boundary  conditions  and 
calibration  and  verification.  This  effort  included  installation  of  five  tide 
gages  at  appropriate  locations  within  the  study  area  (Figure  4).  These  in¬ 
struments  were  maintained  from  23  January  to  12  February  1985  for  a  period  of 
3  weeks  (Table  2).  Tide  elevation  data  were  recorded  at  every  5-min  interval 
during  the  data  collection  period. 

33.  A  harmonic  analysis  of  tide  data  was  conducted  to  remove  data 
spikes  and  high  and  low  frequency  trends,  and  to  subtract  the  mean  from  the 
data  record.  The  mean  water  level  was  referenced  to  NGVD.  Filtered  tide  data 
were  then  analyzed  for  the  amplitude  and  epoch  (phase)  of  their  respective 
tidal  constituents.  This  harmonic  analysis  was  necessary  to  obtain  boundary 
conditions  for  input  to  WIFM.  The  constituents,  amplitudes,  and  epoch 
obtained  from  the  harmonic  analysis  are  presented  in  Table  3-  The  principal 
tidal  constituents  analyzed  were  three  semidiurnal  (M2,  S2,  and  N2)  and  three 
diurnal  (01,  K 1 ,  and  PI).  These  six  major  tidal  constituents  were  selected  to 
provide  boundary  conditions  for  numerical  model  calibration  and  verification. 
Details  of  the  tidal  constituent  analysis  are  reported  by  Schmalz  (1985). 
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Figure  6.  WES  hand-held  current  meter  locations 

31*.  Plates  9-12  show  the  results  of  a  harmonic  analysis  and  the  mea¬ 
sured  water  surface  elevation  records  for  the  offshore,  Rickenbacker  Causeway, 
MacArthur  Causeway  East,  and  Fisher  Island  gages.  These  results  indicate  that 
the  comparison  between  measured  and  calculated  water  elevations  above  NGVD  is 
quite  good,  as  should  as  expected. 

Effect  of  Gulf  Stream  on  Tidal  Circulation 

35.  Coastal  currents  under  the  influence  of  winds,  waves,  and  sea  level 
fluctuations  generate  boundary  currents  such  as  those  of  the  Gulf  Stream. 

Most  of  the  water  in  the  Gulf  Stream  comes  through  Yucatan  Channel  from  the 
Caribbean  Sea  (Stommel  1965).  The  water  is  forced  through  the  long  channel 
caused  by  a  differential  head  of  about  19  cm  (0.62  ft)  between  the  Florida 
Peninsula  on  one  side  and  the  island  of  Cuba  on  the  other.  This  channel  be¬ 
comes  narrower  and  shallower  downstream  and  attains  a  minimum  cross  section 
off  the  Coast  of  Miami  (about  20  miles)  at  longitude  80°30'  and  latitude 
25° 15'  where  the  width  is  about  130  miles  and  the  maximum  depth  is  2,500  ft. 
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The  average  velocity  of  the  Gulf  Stream  in  this  region  is  about  0.3  fps,  and 

Q 

water  mass  transport  is  estimated  to  be  9.2  x  10  cfs  (Stommel  1965).  The 
order  of  magnitude  of  the  maximum  surface  velocity  in  the  channel  off  Miami  is 
about  1.5  fps.  The  axis  of  maximum  surface  velocity  is  not  in  the  center  of 
the  upper  layer  of  the  Gulf  Stream  but  is  displaced  toward  Miami.  This  eccen¬ 
tricity  of  velocity  generates  anticyclonic  vorticity  in  the  upper  layer.  How- 

i 

ever,  this  upper  layer  does  not  extend  into  the  study  area.  (The  study  area 
is  less  than  10  miles  off  the  coast,  while  the  predominant  Gulf  Stream 
influence  is  20  miles  off  the  coast.) 

36.  Sea  level  change  along  the  Miami  coastline  is  minimum  during  July 
and  August  (about  5  cm)  and  maximum  during  October  (8  cm).  The  diurnal  sea 
level  signal  compared  with  the  semidiurnal  signal  is  very  small  at  the  Miami 
coast.  The  magnitude  of  current  (other  than  tide  induced)  generated  from  sea 
level  changes  is  on  the  order  of  7  cm/sec  (0.23  fps)  during  summer  and  about 
5  cm/sec  (0.16  fps)  during  winter  (Blaha  1984).  These  results  indicate  that 
the  contribution  of  sea  level  fluctuations  to  the  observed  east-west  and 
north-south  component  of  current  velocity  is  insignificant. 

37.  The  period  selected  for  numerical  model  calibration  (26  January 
1985)  and  verification  (7  February  1985)  does  not  show  Gulf  Stream  influence. 
Nonexistance  of  Gulf  Stream  Influence  is  further  demonstrated  by  the  close 
agreement  between  the  measured  tides  and  calculated  time  series  of  water  sur¬ 
face  elevations  generated  from  tidal  constituents.  It  is  concluded  that  Gulf 
Stream  currents  have  a  negligible  effect  on  the  overall  tidal  circulation  at 
the  entrance  to  the  port  and  in  the  port  facility  areas  during  most  parts  of 
the  year.  While  it  is  possible  that  the  Gulf  Stream  might  occasionally  mi¬ 
grate  closer  to  the  harbor  entrance,  it  is  impossible  to  predict  these  occur¬ 
rences  or  to  quantify  their  effect  on  tidal  circulation. 
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PART  V:  INPUT  REQUIREMENTS  AND  BOUNDARY  CONDITIONS 


Input  Requirements 


38.  The  numerical  model  WIFM  used  in  this  study  requires  an  input  data 
set  consisting  of  27  card  groups.  Also  associated  with  WIFM  are  preprocessing 
and  postprocessing  codes  which  provide  data  to  or  a  means  of  interpreting  the 
results  of  WIFM.  Care  must  be  taken  in  noting  which  groups  of  data  should  be 
omitted  and  the  number  of  cards  needed  for  these  groups  of  data.  The  com¬ 
plete  input  requirements  for  WIFM  are  explained  in  detail  by  Butler  (in 
preparation) . 


Boundary  Conditions 


39.  A  number  of  boundary  conditions  are  allowed  in  WIFM  which  can  be 
classified  into  three  categories:  open-water  boundaries,  water-land  bound¬ 
aries,  and  subgrid  barrier  boundaries.  A  complete  description  is  given  by 
Butler  (in  preparation). 

Open  boundaries 

40.  The  open  boundary  category  includes  seaward  boundaries  terminating 
the  computational  grid  or  channels  exiting  the  two-dimensional  grid  at  any 
point  in  the  grid  system.  Water  elevations  or  flow  velocities  can  be  pre¬ 
scribed  as  functions  of  location  and  time.  This  information  can  be  input  to 
WIFM  either  as  tabular  data  or  as  constituent  tidal  components  from  which 
water  surface  elevations  can  be  calculated  within  the  code  during  the  time¬ 
marching  process.  In  the  present  study,  tabular  water  levels  were  prescribed 
at  the  open-water  boundary. 

Land/Water  boundaries 


41.  Included  in  the  land/water  boundaries  are  either  fixed  or  variable 
boundaries  to  allow  flooding  and  drying  of  cells.  The  usual  condition  at 
these  boundaries  is  "no-flow"  normal  to  the  boundary  which  is  accomplished  by 
setting  u  =  o  or  v  =  o  at  the  appropriate  cell  face.  Low-lying  terrain 
may  alternately  dry  and  flood  within  a  tidal  cycle  or  surge  history.  Flooding 
in  WIFM  is  simulated  by  making  the  location  of  the  land/water  boundary  a  func¬ 
tion  of  local  water  depth.  Once  the  water  levels  in  adjacent  cells  rise  above 
the  level  of  adjacent  land  height  (possibility  of  flooding),  water  is  allowed 
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to  flow  into  the  "dry"  cell  according  to  a  broad  created  weir  formula  (Reid 
and  Bodine  1968).  When  the  water  level  on  the  dry  cell  exceeds  some  small 
prescribed  value,  the  boundary  face  is  treated  as  open;  and  computations  for 
n  ,  u  ,  and  v  are  made  for  the  cell  (wet  cell).  The  drying  of  cells  is  the 
inverse  process,  and  water  mass  balance  is  conserved  in  these  procedures. 
Subgrid  barriers 

42.  Subgrid  barriers  are  defined  along  cell  faces  and  are  of  three 
types:  exposed,  submerged,  and  overtopping.  Exposed  barriers  permit  no  flow 
across  a  cell  face.  Submerged  barriers  are  simulated  by  controlling  flow 
across  cell  faces  with  the  use  of  a  time-dependent  frictional  coefficient. 
Overtopping  barriers  are  treated  by  using  the  broad-crested  weir  formula  which 
calculates  the  proper  flow  rate  across  the  barrier.  The  barrier's  character¬ 
istics  are  determined  by  its  height  and  water  elevations  in  the  two  adjoining 
cells. 

43.  At  the  open-water  boundaries,  tidal  constituents,  tidal  elevations, 
or  discharge  input  data  may  be  used  to  drive  WIFM.  In  this  study  the  gen¬ 
erated  tidal  constituents  were  specified  at  the  open-water  boundaries.  The 
code  uses  these  constituents  to  generate  time  series  of  surface  elevations  at 
the  boundaries. 


Tidal  Constituents 


Global  grid 

44.  Tidal  constituents  were  specified  as  the  boundary  condition  at  the 
following  global  grid  boundaries:  open  ocean,  Rickenbacker  Causeway,  Mac- 
Arthur  Causeway  East,  MacArthur  Causeway  West,  and  Coast  Guard  station.  Fig¬ 
ure  1  shows  these  boundaries. 

Refined  grid 

45.  Tidal  constituents  were  specified  as  the  boundary  condition  at  the 
following  refined  grid  boundaries:  MacArthur  Causeway  East,  MacArthur  Cause¬ 
way  West,  and  Coast  Guard  station.  Hydrodynamic  values  were  transferred  from 
the  global  grid  to  the  refined  grid  at  the  open-ocean  boundary,  northern 
boundary,  and  at  the  southern  boundary,  including  Rickenbacker  Causeway. 
Hydrodynamic  data  transfer  was  accomplished  by  linking  computational  grids  of 
the  global  grid  with  the  refined  grid  as  described  in  Part  III  of  this  report. 
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PART  VI:  NUMERICAL  MODEL  CALIBRATION  AND  VERIFICATION 


Global  Grid  Calibration 

46.  WIFM  was  calibrated  against  the  time  series  of  water  surface  ele¬ 
vations  developed  from  the  tidal  constituents  generated  from  the  harmonic 
analysis  of  prototype  data  for  water  surface  elevations  at  Rickenbacker  Cause¬ 
way,  MacArthur  Causeway  East,  Fisher  Island,  and  an  offshore  gage.  Six  pre¬ 
dominant  tidal  constituents  (M2,  S2,  N2,  K^,  0^,  and  1^)  were  used  for  the 
numerical  model  calibration.  These  constituents  were  selected  because  they 
provided  a  good  representation  of  tide  for  the  Miami  coast. 

47.  The  capability  of  the  numerical  model  to  accurately  predict  pro¬ 
totype  data  depends  on  many  factors:  (a)  boundary  specified  conditions, 

(b)  cell  size  in  the  grid,  (c)  water  depths  assigned  to  each  cell,  and 

(d)  Manning's  roughness  coefficient  for  each  cell.  A  detailed  description 
of  the  sensitivity  of  WIFM  results  to  these  parameters  can  be  found  in 
Leenknecht  et  al.  (1984).  In  this  study  the  elevations  assigned  to  each  cell 
and  the  Manning's  coefficient  were  adjusted  to  match  predicted  tides  with 
prototype  data. 

48.  The  WIFM  calibration  period  was  chosen  to  be  in  the  neap  tidal 
period  beginning  at  midnight  26  January  1985.  This  period  was  selected 
because  the  agreement  between  measured  and  calculated  tides  is  the  best,  and 
accurate  velocity  measurements  were  available.  In  addition,  a  majority  of  the 
field  stations  have  recorded  prototype  data  during  this  period  (Figure  5). 
Tidal  constituents  (magnitude  and  phase)  were  specified  at  Rickenbacker  Cause¬ 
way,  MacArthur  Causeway  East,  Fisher  Island,  Miami  Marina,  and  the  offshore 
boundary  of  the  computational  grid  (see  Figure  1).  Tidal  constituents  for 
Miami  Marina  were  obtained  from  NOAA.  The  tidal  constituents  specified  at  the 
different  gage  locations  within  the  global  grid  served  as  boundary  conditions 
for  the  numerical  model  calibration. 

49.  A  60-sec  time-step  was  used  in  WIFM,  and  26  hr  of  the  tide  were 
simulated  at  20  locations  within  the  global  grid.  Four  stations  (Rickenbacker 
Causeway,  MacArthur  Causeway,  Fisher  Island,  and  offshore)  were  selected  for 
comparisons  with  WIFM  computations  during  the  entire  calibration  period.  The 
number  of  cells  and  water  depth  for  cells  were  adjusted  to  adequately  describe 
flow  conditions  at  the  internal  boundaries  (MacArthur  Causeway  East  and  West, 


Coast  Guard  station,  and  Rickenbacker  Causeway). 

50.  Plates  13-16  compare  predicted  and  measured  tides  at  the  bound¬ 
aries.  Plate  17  presents  a  multiplot  of  calculated  tidal  elevations  for 
boundary  gages.  Plate  18  compares  measured  tide  at  a  nonboundary  gage  (Fisher 
Island)  and  calculated  results.  The  model  reproduced  the  offshore  tide  with 
reasonable  accuracy.  Model  performance  at  other  locations  is  good  with  ampli¬ 
tude  variation  less  than  0.1  ft  and  a  slight  phase  shifting  caused  by  channel 
and  shoaling  features  of  the  study  area. 

51.  Plates  19  and  20  present  the  results  of  WIFM  compared  to  the  pro¬ 
totype  current  measured  by  NOAA.  It  is  seen  that  velocity  records  of  peak 
flood  and  ebb  flows  and  calculated  peak  flood  and  ebb  flows  match  quite  well 
for  the  stations  shown.  Because  of  the  failure  of  several  NOAA  field  sta¬ 
tions,  only  selective  comparisons  are  shown. 

52.  Plates  21  and  22  show  computed  current  fields  over  a  portion  of  the 
grid  in  the  form  of  vector  plots  for  a  flood  and  an  ebb  flow,  respectively. 
Simulation  time  and  vector  scale  are  given  at  the  bottom  of  each  plate.  Tidal 
circulation  patterns  appear  quite  reasonable  in  direction  and  in  magnitude. 

53.  The  above  results  indicate  that  WIFM's  ability  to  describe  the- 
hydrodynamics  of  the  area  of  interest  is  quite  good.  The  model  quite 
accurately  reproduced  the  observed  tide  and  velocity  at  all  stations  where 
measurements  were  available. 

Global  Grid  Verification 

54.  Verification  of  WIFM  consists  of  its  ability  to  reproduce  accu¬ 
rately  the  hydrodynamics  of  the  area  of  interest  without  adjusting  model 
parameters.  To  test  this  accuracy,  a  period  on  7  February  1985  was  simulated 
by  WIFM.  This  verification  period  was  selected  because  extreme  spring  tidal 
events  were  recorded  along  the  Miami  coast,  and  they  have  significant  effects 
on  tidal  circulation  in  the  study  area.  In  addition,  a  majority  of  the  field 
stations  have  accurate  prototype  velocity  data  measured  during  this  period 
(Figure  5) . 

55.  To  verify  the  model,  computation  in  WIFM  began  at  midnight  on 

7  February  1985  using  the  predominant  tidal  constituents  described  earlier. 

The  simulation  continued  for  a  period  of  26  hr.  A  60-sec  time-step  was  used 
in  the  model. 
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56.  Computed  tidal  elevation  and  flow  velocities  were  compared  with 
prototype  data  at  various  locations  over  the  grid.  Plates  23-26  compare  cal¬ 
culated  and  measured  tidal  elevations  at  the  boundary  gages.  Plate  27  shows  a 
multiple  plot  of  WIFM-calculated  tidal  elevations  of  boundary  gages.  Plate  28 
compares  the  measured  and  calculated  tides  at  a  nonboundary  gage  (Fisher 
Island) . 

57.  Plates  29-31  show  the  measured  and  computed  flow  velocities  which 
were  obtained  for  the  model  verification  period.  Comparison  between  the  mea¬ 
sured  and  calculated  currents  is  quite  good  except  at  NOAA  Stations  10  and  17. 

58.  Figure  7  was  constructed  to  show  a  comparison  of  WIFM  current 
computations,  prototype  data  measured  by  NOAA,  and  hand-held  data  collected  by 
WES  at  NOAA  Station  10.  The  comparison  between  WIFM  computation  and  the  hand¬ 
held  data  is  quite  good,  while  prototype  data  from  the  NOAA  gage  show  some 
differences.  Plate  32  displays  a  good  comparison  of  WIFM  computations  with 
hand-held  data  at  various  locations  along  the  grid  (Figure  2).  NOAA's  pro¬ 
totype  data  are  not  available  at  these  locations. 

59.  Plates  33  and  3*1  show  computed  current  fields  in  the  form  of  vector 
plots  for  a  selected  portion  of  the  grid.  These  plots  represent  peak  flood- 
and  ebb-tidal  flow  through  the  entrance.  Simulation  time  and  vector  scale  are 
given  at  the  bottom  of  each  plate. 

Refined  Grid  Verification 

60.  Verification  of  the  refined  numerical  model  was  accomplished  by 
using  the  refined  grid  boundary  conditions.  The  minimum  cell  dimension  was 
restricted  to  120  ft.  A  60-sec  time-step  gave  simulation  stability  and  pro¬ 
vided  sufficient  resolution  to  accurately  describe  tidal  circulation  in  the 
area  of  interest.  To  keep  the  expenses  of  verification  to  a  minimum,  WIFM 
simulated  16  hr  of  tide  during  spring  tide  (7  February  1985). 

61.  Plates  35-38  compare  calculated  and  measured  tidal  elevations  of 
the  boundary  gages.  These  results  show  insignificant  variation  in  amplitude 
and  phase  between  the  measured  and  calculated  water  surface  elevations. 

Plate  39  compares  the  measured  tide  at  a  nonboundary  gage  (Fisher  Island)  and 
WIFM-calculated  results.  Peak  amplitude  variation  between  measured  and  cal¬ 
culated  results  is  less  than  0.15  ft  with  little  phase  shifting. 

62.  Plates  40  and  41  compare  WIFM-calculated  current  magnitudes  and 
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directions  with  prototype  data  measured  by  NOAA.  The  model  reproduced  pro¬ 
totype  velocity  records  with  proper  phasing  of  the  peak  flood  and  ebb  flows. 

63.  Based  on  the  results  described  above,  it  was  concluded  that  the 
refined  numerical  model  accurately  simulated  the  hydrodynamics  of  the  area  of 
interest,  and  calibration  of  the  refined  numerical  model  was  not  required. 
Therefore,  no  hydrodynamic  simulations  were  performed  for  neap  tide  conditions 
(25  January  1985). 


. . 


PART  VII:  EVALUATION  OF  ENTRANCE  CHANNEL  CROSS-CURRENT 
PROBLEMS  AND  INNER-HARBOR  CHANNEL  NAVIGATION  PROBLEMS 


Cross-Current  Problems 


64.  Production  tests  were  conducted  to  assist  CESAJ  in  assessing  the 
cross-current  problems  at  the  entrance  to  the  Port  of  Miami.  The  investiga¬ 
tion  involved  a  numerical  modeling  effort  to  identify  the  problems  and  to 
evaluate  the  impact  of  proposed  corrective  alternatives  (structural  and  non- 
structural)  on  the  hydrodynamic  regime  near  the  port  entrance.  In  addition 
to  the  entrance  problem,  CESAJ  is  interested  in  evaluating  the  hydrodynamic 
impact  of  proposed  changes  (widening  and  deepening  of  harbor  channels)  to  port 
facilities.  Therefore,  the  numerical  model  test  results  cover  a  substantial 
portion  of  the  Port  of  Miami  and  adjacent  areas.  On  the  basis  of  the  results 
presented  for  refined  numerical  model  verification,  it  is  concluded  that  WIFM 
correctly  simulates  Port  of  Miami  hydrodynamics  and  that  the  refined  numerical 
model  can  therefore  be  used  for  evaluating  cross-current  problems  at  the  en¬ 
trance  to  the  port  and  navigation  problems  in  the  inner  harbor  channels. 

No  change  in  topography 

65.  For  the  the  no  change  in  topography  test,  and  subsequent  tests,  the 
testing  period  for  hydrodynamic  scenarios  was  7  February  1985  (spring  tide 
event),  the  refined  numerical  grid  was  used,  and  the  forcing  boundary 
conditions  were  identical  to  those  used  in  the  refined  grid  verification.  The 
above  testing  period  was  selected  because  extreme  spring  tidal  events,  which 
have  significant  effects  on  tidal  circulation  in  the  study  area,  were  recorded 
along  the  Miami  coast. 

66.  A  numerical  simulation  was  conducted  with  WIFM  to  determine  the 
coastal  current  regime  and  harbor  hydrodynamics  for  the  existing  topography. 
The  simulation  allowed  comparison  of  subsequent  test  results  obtained  for 
various  structural  and  nonstructural  modification  plans.  Hydrodynamic  results 
were  simulated  for  one  tidal  cycle  P2.8  hr)  using  a  60-sec  time-step.  This 
testing  period  included  flood  and  ebb  flows  through  the  entrance  channel. 
Plates  42  and  43  show  computed  current  field  for  the  flood  and  ebb  flows  in 
the  form  of  vector  plots  for  a  selected  portion  of  the  grid.  Simulation  time 
and  vector  scale  are  given  at  the  bottom  of  each  plate. 
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Structural  changes 

67.  Modified  Jetty  configurations  were  simulated  using  existing  ba¬ 
thymetry  features.  The  structural  changes  tested  included  four  alternative 
Jetty  modification  plans  to  alleviate  problems  with  strong  cross  currents  at 
the  ocean  entrance.  These  alternative  plans  are  given  in  Table  4  (Cases  2-5). 

68.  Plates  44  and  45  show  computed  current  field  for  the  flood  and  ebb 
flows  in  the  form  of  vector  plots  for  the  second  test  (Table  4).  A  compari¬ 
son  between  test  results  obtained  for  the  existing  topography  (Plates  42  and 
43)  and  the  second  test  (Plates  44  and  45)  indicates  that  the  strong  cross 
currents  generated  at  the  tip  of  the  south  Jetty  due  to  unequal  Jetty  lengths 
were  eliminated  when  the  jetties  were  of  equal  length.  This  finding  clearly 
indicates  that  the  effects  of  cross  currents  on  navigation  are  greater  when 
the  jetties  are  of  unequal  lengths.  A  possible  solution  for  alleviating  cross 
currents  at  the  entrance  is  even  length  jetties. 

69.  Plates  46  and  47  present  computed  current  field  for  the  flood  and 
ebb  flows  which  were  obtained  for  the  third  test  (Table  4).  In  this  case  the 
jetties  were  of  even  length,  and  the  north  jetty  was  flared.  The  purpose  of 
providing  flare  at  the  north  jetty  was  to  demonstrate  the  effect  of  confined 
flow  on  tidal  circulation  at  the  entrance.  A  comparison  between  Plates  42 
and  43  and  46  and  47  shows  that  the  magnitude  of  currents  near  the  tip  of  the 
north  Jetty  was  decreased  because  of  the  addition  of  flare  at  the  north  jetty. 

70.  Results  for  the  fourth  test  (Table  4)  are  shown  in  Plates  48  and  49 
for  the  flood  and  ebb  flows,  respectively.  For  this  test  the  Jetties  were  of 
even  length,  and  the  south  jetty  was  flared.  Results  showed  that  providing 
flare  at  the  tip  of  the  south  Jetty  reduces  the  magnitude  of  current  at  the 
flared  end.  However,  strong  cross  currents  at  the  tip  of  the  unflared  jetty 
were  not  evident. 

71.  Plates  50  and  51  present  computed  current  field  for  the  fifth  test 
(Table  4).  In  this  case,  Jetties  were  of  even  length  and  flared  at  the  jetty 
ends.  As  anticipated,  flare  at  the  Jetty  ends  reduced  the  magnitude  of  cur¬ 
rents  and  alleviated  cross  currents.  These  wider  Jetty  flares  produced 
smaller  velocities  at  the  vicinity  of  the  flares  and  negligible  velocity 
changes  along  Government  Cut.  This  combination  of  even  jetty  length  and 
flare  at  Jetty  ends  appears  to  be  a  corrective  measure  for  alleviating  the 
problem  with  strong  cross  currents  at  the  ocean  entrance  to  the  port  of  Miami. 


m 


vn 


■  *•>  1 


X\sv- 


31 


* , » ,  -  ***«• .  * ,  *.  **.  • ,  *■.  •*  /  /  .* 


r.^  V.  »*-  < 


■v 

.  j T-  »  ^  -  m 


VX- -/-.v 

I  ■ 


vVvl 

\*Vv 

V,  v, 

c  vv,’ 

■m 


.v:v>:-:vX-h 


Nonstructural  changes 

72.  In  the  event  jetty  modification  is  precluded  because  of  lack  of 
funding,  consideration  should  be  given  to  nonstructural  modifications  (widen¬ 
ing  and  deepening  of  navigation  channels).  This  alternative  may  be  less 
costly  than  structural  changes  while  providing  the  additional  room  required 
for  safe  maneuvering  of  large  vessels  through  strong  cross  currents  at  the 
ocean  entrance. 

73.  Test  6  (Table  4)  was  conducted  with  WIFM  to  examine  tidal  circula¬ 
tion  changes  resulting  from  a  hole  at  the  tip  of  the  north  jetty.  The  size  of 
the  hole  and  the  water  depth  in  the  hole  are  shown  in  Plate  52.  Plates  52  and 
53  show  computed  current  fields  for  the  flood  and  ebb  conditions,  respec¬ 
tively.  A  comparison  of  Plate  42  with  Plates  52  and  53  shows  that  the  magni¬ 
tude  of  current  in  the  hole  is  reduced  because  of  increased  water  depth  in  the 
hole. 


Navigation  Problems 

74.  CESAJ  provided  a  list  of  projected  water  depths  in  the  harbor  chan¬ 
nels  along  with  information  that  defines  the  proposed  channel  and  turning 
basin  widths  to  be  tested  using  WIFM.  Table  5  lists  the  maximum  project  water 
depths  which  were  used  with  WIFM  for  channel  deepening.  Figure  8  shows  as¬ 
signed  water  depths  (which  are  circled)  at  various  locations  in  the  harbor  and 
the  boundary  lines  of  widened  channels. 

Nonstructural  changes 

75.  Test  7  (Table  4)  was  conducted  with  WIFM  to  examine  the  nonstruc¬ 
tural  impact  on  tidal  circulation  at  the  entrance  and  in  the  harbor  channels. 
The  nonstructural  changes  included  channel  deepening  and  widening  as  shown  in 
Figure  7.  Table  5  lists  the  maximum  project  water  depths  used  for  channel 
deepening.  All  tests  were  conducted  for  the  spring  tide  events  (7  February 
1985).  Plates  54  and  55  present  computed  tidal  currents  for  the  flood  and  ebb 
flows,  respectively.  A  comparison  of  test  results  obtained  for  the  existing 
topography  (Plates  42  and  43)  and  the  results  (Plates  54  and  55)  of  this  test 
indicate  that  the  overall  tidal  circulation  in  the  harbor  does  not  change  sig¬ 
nificantly,  and  the  magnitude  of  the  current  increases  slightly  in  areas  sub¬ 
jected  to  channel  deepening  and  widening. 


Boundary  lines  of  proposed  channel  deepening  and  widening 


Structural  and  nonstructural  changes 


76.  For  practical  purposes,  sometimes  a  combination  of  structural  and 
nonstructural  changes  is  warranted  for  a  coastal  project.  Table  4  lists  two 
of  these  modification  plans  (Tests  8  and  9)  tested  in  WIFM. 

77.  In  Test  8,  the  Jetty  lengths  were  made  even,  and  the  nonstructural 
changes  included  deepening  at  the  jetty  ends,  as  shown  in  Plate  56.  No  other 
changes  were  made  in  the  harbor  channels  (existing  topography  was  used). 

Plates  56  and  57  present  computed  current  field  for  the  flood-  and  ebb-tidal 
events,  respectively.  A  comparison  of  this  test  to  test  results  obtained  for 
the  second  test  indicates  that  the  strong  cross  currents  at  the  entrance  to 
the  port  were  reduced  by  deepening  at  the  inlet  throat.  Results  in  this  com¬ 
parison  also  indicate  the  holes  have  no  impact  on  currents  in  the  main  chan¬ 
nel.  Consequently,  widening  and  deepening  the  entrance  channel  appear  to  be 
appropriate  to  mitigate  cross-current  problems  at  the  entrance. 

78.  The  final  test  (Test  9)  conducted  in  WIFM  included  even  length 
jetties,  flared  jetty  ends,  and  deepening  and  widening  of  the  inlet  and  harbor 
channels  (similar  to  Test  7).  Plates  58  and  59  present  the  computed  tidal 
circulation  structure  at  the  entrance  and  port  facility  areas  for  the  flood- 
and  ebb-tidal  flows,  respectively.  It  is  seen  that  the  overall  circulation  in 
the  harbor  did  not  change  (compare  results  of  Tests  7  and  9)  with  the  struc¬ 
tural  changes  provided  at  the  entrance  to  the  port.  However,  the  problem  with 
strong  cross  currents  at  the  entrance  to  the  Port  of  Miami  was  alleviated. 

79.  Table  6  shows  the  relative  changes  in  peak  flood  and  ebb  velocities 
at  selected  points  with  and  without  harbor  improvements.  An  examination  of 
results  presented  in  Table  6  indicates  that  harbor  improvements  would  not  have 
significant  impact  on  current  flow  through  the  harbor. 


PART  VIII:  PRODUCTION  RUNS  FOR  INPUT  TO  VESSEL  SIMULATION  STUDY 


80.  CAORF  is  conducting  a  vessel  simulation  study  sponsored  by  CESAJ 
for  the  Port  of  Miami.  This  study  requires  the  refined  numerical  model  test 
results  obtained  for  the  existing  navigation  channels  with  existing  jetties, 
with  Jetties  extended  to  be  equal  and  flared,  and  for  deepened  and  widened 
channels.  These  various  tests  are  listed  in  Table  7.  For  each  test,  the 
input  data  into  the  vessel  simulation  study  were  obtained  during  two  different 
peak  tidal  flood  conditions  (at  hour  4),  one  in  the  absence  of  wind  and  one 
with  a  constant  wind  of  20  knots  from  the  northeast. 

Global  Model  Test  Runs 

81.  Two  test  runs  with  the  global  numerical  model  were  conducted  in 
WIFM  using  the  global  grid.  For  these  tests,  and  subsequent  ones  (global  or 
refined),  the  testing  period  for  hydrodynamic  scenarios  was  7  February  1985 
(at  peak  spring  tide),  and  the  forcing  boundary  conditions  without  the  wind 
were  identical  to  those  used  in  the  global/refined  grid  verification.  The 
first  test  dealt  with  the  existing  topography  and  the  boundary  conditions  de¬ 
scribed  above.  The  second  test  was  similar  to  the  first  but  included  the 
appropriate  wind. 

82.  Plate  60  shows  computed  current  field  vector  plots  of  the  first 
test  for  the  maximum  flood  flows  for  a  selected  portion  of  the  grid.  Simula¬ 
tion  time  and  vector  scale  are  given  at  the  bottom  of  the  plate.  The  hydro- 
dynamic  results  from  the  global  model  were  saved  for  input  into  the  fine  reso¬ 
lution  model. 

83.  The  second  test  included  a  wind  of  20  knots  from  the  northeast, 
which  was  superimposed  on  the  existing  topography  of  the  study  area  and  the 
previously  used  boundary  conditions  for  the  first  test.  Plate  61  presents 
the  computed  current  vector  fields  for  the  peak  flood  flows.  A  comparison  of 
Plates  60  and  61  indicates  that  the  effect  of  storm  winds  at  the  entrance  to 
the  Port  of  Miami  is  significant.  These  plots  show  a  distinct  difference  in 
the  magnitude  and  direction  of  the  velocity  vectors  at  the  entrance  to  the 
Port.  The  hydrodynamic  results  for  the  second  test  were  saved  for  input  into 
the  refined  model. 


Refined  Model  Test  Runs 

84.  A  numerical  test  (Test  3,  Table  7)  was  conducted  in  the  refined 
numerical  model  to  determine  the  coastal  current  regime  and  harbor  hydro¬ 
dynamics  for  the  existing  topography.  The  forcing  boundary  conditions  at 
MacArthur  Causeway  East,  MacArthur  Causeway  West,  and  the  Coast  Guard  station 
were  identical  to  those  used  in  the  refined  model  verification.  However,  the 
boundary  conditions  specified  at  Rickenbacker  Causeway  and  at  the  offshore 
boundary  of  the  computational  grid  were  the  hydrodynamic  results  obtained  from 
the  global  grid  for  peak  flood  flow  through  the  entrance  channel.  Hydro- 
dynamic  scenarios  were  obtained  for  hour  4  using  a  60-sec  time-step.  Plate  62 
shows  computed  tidal  currents  for  the  flood  flow  through  the  entrance  channel 
in  the  form  of  vector  plots  for  a  selected  portion  of  the  grid.  Simulation 
time  and  vector  scale  are  given  at  the  bottom  of  the  plate. 

85.  Results  for  the  fourth  test  (Table  7)  are  shown  in  Plate  63  for  the 
peak  flood  flow  conditions.  For  this  test,  the  Jetties  were  extended  to  be 
equal  in  length,  and  the  boundary  conditions  used  in  WIFM  were  for  calm  con¬ 
ditions. 

86.  Plate  64  presents  computed  current  field  for  the  peak  flood  flow 
through  the  entrance  for  the  fifth  test  (Table  7).  In  this  case,  the  Jetties 
are  set  equal  in  length  and  flared,  and  the  entrance  and  harbor  channels  are 
widened  and  deepened.  The  effect  by  winds  on  the  overall  tidal  circulation  of 
the  study  area  was  not  included  in  this  test. 

87.  Plates  65,  66,  and  67  present  corresponding  plots  which  were  ob¬ 
tained  for  a  constant  wind  20  knots  from  the  northeast.  A  comparison  of  these 
results  (Plates  65-67)  to  test  results  (Plate  62-64)  obtained  for  the  3rd, 

4th,  and  5th  tests  indicates  that  the  magnitude  and  direction  of  currents  at 
the  entrance  to  the  port  were  intensified  in  the  presence  of  the  20-knot  wind 
from  the  northeast.  Table  8  lists  the  magnitude  and  direction  (measured  from 
true  north)  of  velocity  (with  and  without  wind  effects)  at  selective  numerical 
grid  point  locations.  The  model  indicates  that  the  magnitude  of  the  velocity 
at  the  tip  of  the  north  jetty  was  increased  by  about  0.5  fps,  and  the  direc¬ 
tion  changed  by  about  10  deg  because  of  the  wind.  This  effect  is  seen  for 
Cases  2  and  4  (Table  8).  For  Case  6,  the  magnitude  of  velocity  was  increased 
by  0.2  fps,  and  the  change  in  direction  was  insignificant.  At  the  south 
Jetty,  the  magnitude  and  direction  of  the  velocity  (Table  8,  Cases  1  and  2) 
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were  changed.  However,  when  the  Jetties  were  extended  to  equal  length  (Table 
8,  Cases  3  and  4),  the  magnitude  of  velocity  was  increased  by  about  0.2  fps 
and  the  direction  was  unchanged.  Results  for  Cases  5  and  6  indicate  that  the 
magnitude  of  velocity  increased  by  about  0.2  fps  at  the  three  locations 
selected  for  model  analysis,  while  the  direction  of  velocity  remained  prac¬ 
tically  the  same.  The  magnitude  of  velocity  at  midchannel  increased  by  about 
0.2  fps  for  all  cases.  For  Cases  1-4,  the  direction  of  velocity  changed  by 
about  5  deg  because  of  the  wind.  However,  no  change  in  direction  was  evident 
for  Cases  5  and  6. 


PART  IX:  SUMMARY  AND  CONCLUSIONS 


88.  WES  applied  WIFM  in  a  two-phase  numerical  modeling  study  to  evalu¬ 
ate  tidal  current  regime  at  the  entrance  to  the  Port  of  Miami  and  to  examine 
the  impact  of  the  proposed  port  facility  changes  on  tidal  circulation  in  the 
harbor.  The  first  phase  consisted  of  a  global  model  and  the  second  phase  a 
fine  resolution  model.  The  global  model  was  calibrated  against  tidal  con¬ 
stituents  generated  from  a  harmonic  analysis  of  the  prototype  data.  WIFM- 
predicted  neap  tide3  (26  January  1985)  were  compared  against  the  measured 
data.  The  comparison  was  good. 

89.  The  global  and  refined  numerical  models  were  verified  during  spring 
tides  (7  February  1985).  Computed  water  surface  elevations  above  NGVD  and 
flow  velocities  were  compared  with  the  field  data  at  several  locations  across 
the  grid.  These  comparisons  were  good.  It  was  concluded  that  the  refined 
numerical  model  can  accurately  describe  the  hydrodynamics  of  the  area  of  in¬ 
terest  and  that  calibration  of  this  model  is  not  required. 

90.  Nine  alternative  plans  (structural,  nonstructural ,  and  the  combi¬ 
nation  of  both)  were  tested  in  the  refined  numerical  model.  The  testing  . 
period  for  all  cases  was  7  February  1985  (spring  tide). 

91.  Structural  alternatives  at  the  entrance  to  the  port  were  found  to 
effectively  reduce  cross  currents  when  the  jetties  were  made  equal  length.  In 
addition,  flared  Jetties  of  equal  length  provided  more  effective  results. 

92.  The  nonstructural  modification  included  channel  deepening  and  wid¬ 
ening  at  the  entrance  to  the  port  and  in  the  port  facility  areas.  Hydro- 
dynamic  scenarios  obtained  for  this  case  indicated  that  the  overall  tidal 
circulation  in  the  harbor  did  not  change  significantly  and  tnat  the  magnitude 
of  current  increased  slightly  in  areas  subjected  to  channel  deepening  and 
widening. 

93.  The  structural  and  nonstructural  changes  included  Jetty  extensions 
as  well  as  deepening  and  widening  inlet  and  harbor  channels.  The  hydrodynamic 
results  obtained  for  Case  9  (Table  4)  (flared  jetties  of  equal  length  and 
changes  made  to  inlet  and  harbor  channels)  were  effective  in  alleviating 
strong  cross  currents  at  the  entrance. 

94.  The  numerical  modeling  effort  required  for  the  Port  of  Miami  Ship 
Simulation  Study  was  investigated  during  two  different  peak  tidal  flood 
conditions  (spring  tide  at  hour  4),  one  in  the  absence  of  wind  and  one  with  a 
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constant  wind  of  20  knots  from  the  northeast.  A  total  of  eight  test  runs  (two 
global  and  six  refined)  was  conducted  with  WIFM  (Table  7). 

95.  Hydrodynamic  results  obtained  for  these  tests  were  provided  for 
input  information  into  the  vessel  simulation  study. 
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Current  Meter 

Current 

Meter 

Identi- 

Data  Collection 

Collec- 
t  ion 
Period 

Position 

f ication 

Latitude 

Longitude 

Began 

Ended 

days 

Nearshore,  north  of 
Government  Cut 

77.CM1 

25°46.2' 

80°7.7' 

1-23-85 

2-12-85 

21 

Near  tip  of  north 
Jetty 

76, CM2 

25°45.8' 

80°7.4’ 

1-23-85 

2-12-85 

21 

Near  tip  of  south 
Jetty 

80, CM3 

25°45.4' 

80o7.6’ 

1-27-85 

2-12-85 

17 

Norris  Cut 

79,CM4 

25°45.5' 

80°8.7' 

1-23-85 

2-12-85 

17 

Near  Outer  Bar 

78,G01 

25°45.7' 

80°5.9’ 

1-23-85 

2-25-85 

34 

Table  2 

Summary  of  Prototype  Tide  Data,  WES  Effort 


Tide  Gage 
Location 

Tide 

Gage 

Identi¬ 

fication 

Latitude 

Longitude 

Data  Collection 
Began  Ended 

Collection 
Period , 
days 

Fisher  Island 

TG  49 

25°45.3' 

80°8.5' 

1-22-85 

2-12-85 

22 

Rickenbacker 

TG  52 

25° 14.7' 

80° i i . r 

1-22-85 

2-12-85 

22 

Causeway 

Offshore 

TG  72 

25°45.8' 

80°8.0' 

1-22-85 

2-12-85 

22 

MacArthur 

TG  85 

25°47.3’ 

80° 10.9.5' 

1-22-85 

2-12-85 

22 

Causeway ,  W 

MacArthur 

TG  108 

25°46.4' 

80°8.7' 

1-22-85 

2-07-85 

17 

Causeway ,  E 
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Table  3 

Summary  of  Tidal  Constituents  Generated  from  Harmonic  Analysis 


Gage  Location 


Constituents 


Amplitude 

ft 


Epoch 
de 


Root-Mean 
Square  Error 


Fisher  Island 

M2 

1 .04797 

226.16911 

S2 

0.10315 

310.35679 

N2 

0.22081 

217.81348 

0.20096 

K1 

0.14593 

177.94606 

01 

0.09358 

188.681 12 

K2 

0.11529 

153.59656 

Rickenbacker 

Causeway 

M2 

0.92165 

240.21394 

S2 

0.10950 

305.95851 

N2 

0.16834 

236.34916 

0.219685 

K1 

0.12809 

188.70476 

01 

0.09437 

199.32917 

K2 

0.07310 

189.35635 

Offshore 

M2 

1.18837 

216.75510 

S2 

0.06828 

159.38349 

N2 

0 . 28860 

204.94294 

0.202639 

K1 

0.15418 

173.94872 

01 

0. 10946 

180. 10808 

K2 

0.22850 

174.89292 

MacArthur 

Causeway,  E 

M2 

1.01556 

233.43168 

S2 

0.20003 

291 . 18362 

N2 

0.21499 

230.57715 

K 1 

0. 13321 

184.66181 

01 

0.09421 

189.17045 

0.175297 

K2 

0.00889 

203.56246 

* 

Miami  Marina 

M2 

1 .00000 

245.27000 

S2 

0. 17300 

227.25000 

N2 

0.20000 

277.85000 

K 1 

0.10700 

185.99000 

01 

0.08900 

209.09000 

K2 

0.05000 

273.39000 
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Table  4 

Port  of  Miami  Numerical  Model  Study,  Refined  Grid 
Production  Test  Cases 
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Modification  Plan 


Tests 


None 

Structural 

Structural 

Structural 

Structural 


Nonstructural 


Nonstructural 


Structural  and 
nonstructural 


Structural  and 
nonstructural 


Existing  topography 

North  and  south  jetties  of 
equal  length 

North  and  south  Jetties  of 
equal  length  and  flared 
at  the  north  Jetty 

North  and  south  Jetties  of 
equal  length  and  flared 
at  the  south  jetty 

North  and  south  jetties  of 
equal  length  and  flared 
at  Jetty  ends 

Existing  topography  and  a 
hole  provided  at  the  tip  of 
the  north  jetty 

Deepening  and  widening  of 
harbor  channels  and  turning 
basins,  ocean  entrance  channel 
and  outer  bar  channel 

Existing  topography,  both 
jetties  of  equal  length, 
throat  width  increased  at 
the  entrance  on  either  side 
of  Jetty,  and  water  depth 
increased  to  47  ft 

Both  jetties  extended  to  same 
length,  flared  at  jetty  ends, 
and  deepening  and  widening 
of  harbor  channels,  ocean 
entrance  channel  and  outer 
bar  channel 
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Table  5 


Maximum  Project  Water  Depths  Used  for  Channel  Deepeninc 
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Global 

Global 

Refined 

Refined 

Refined 


Refined 

Refined 

Refined 


None 

None 

None 

Structural 


Structural  and 
nonstructural 


None 

Structural 


Structural  and 
nonstructural 


Existing  topography  without  wind 

Existing  topography  with  wind 

Existing  topography  without  wind 

Existing  topography,  even  length 

Jetties  without  wind 

Even  length  jetties,  Jetties 
flared  at  ends,  and  deepen¬ 
ing  and  widening  of  harbor 
channels  without  wind 

Existing  topography  with  wind 

Existing  topography,  even  length 
Jetties,  with  wind 

Even  length  jetties,  jetties 
flared  at  ends,  and  deepening 
and  widening  of  harbor  chan¬ 
nels,  with  winds 


*  For  peak  flood  conditions  with  and  without  20-knot  northeasterly  winds. 
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APPENDIX  A:  NOTATION 


a,b,c 

A,B 

C 

d 

f 

g 

h 

k 

n 

R 


t 

u 

U 

v 

V 

x,y 

x 

z 

At 
Ax,  Ay 
e 

n 


T 
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Regional  constants  derived  from  stretching  transformation  of 
coordinate  system 

Coefficient  matrices 

Chezy  friction  coefficient 

Total  depth  of  water  column  d  =  n  -  h 

Coriolis  parameter 

Acceleration  due  to  gravity 

Local  ground  (cell)  elevation  above  datum 

Integer  time-step  counter 

Dimensionless  parameter  used  to  characterize  stability  criterion 

Rate  of  water  volume  change  in  the  system  (for  example  through 
rainfall  or  evaporation) 

Time 

Vertically  averaged  water  velocity  in  x-direction 

Matrix  consisting  of  n,  u,  and  v,  as  functions  of  x,  y,  and  t 

Vertically  averaged  water  velocity  in  x-direction 

Largest  velocity  encountered  at  a  computational  cell 

Cartesian  coordinate  system  axes 

Smaller  value  of  x  and  y 

Computational  grid  lines  defined  by  positive  integer  value 
Time-step 

Length  of  computational  cell  in  x-  and  y-direction 

Eddy  viscosity  coefficient 

Water  surface  elevation  above  datum 

Hydrostatic  water  elevation  due  to  atmospheric  pressure 
differences 

Two-dimensional  differences  operators 
Air  density 

Surface  stress  of  wind 
Intermediate  time-step  level 

Positive  integer  representing  computational  grid  line 
Centered  difference  operators 
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